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ON THE DYNAMICS OF A THIRD ORDER NEWTON’S 
APPROXIMATION METHOD 

AURELIAN GHEONDEA AND MEHMET EMRE §AMCI 


Abstract. We provide an answer to a question raised by S. Amat, S. Busquier, S. 
Plaza on the qualitative analysis of the dynamics of the third order Newton type 
approximation function Mf, by proving that for functions / twice continuously dif¬ 
ferentiable and such that both / and its derivative do not have multiple roots, with 
at least four roots and infinite limits of opposite signs at ±oo, Mf has periodic points 
of any prime period and that the set of points a at which the approximation sequence 
(My (a))„ 6 N does not converge is uncountable. In addition, we observe that in their 
Scaling Theorem analyticity can be replaced with differentiability. 


1. Introduction 

The classical Newton’s Approximation Function Nf(x ) = x — f(x)/f'(x), for nu¬ 
merical approximation of roots of (nonlinear) functions /, under certain conditions of 
smoothness and distribution of roots and critical points, has a second order speed of 
convergence and, until now, it is considered as one of the most useful and reliable iter¬ 
ative method of this kind. S. Amat, S. Busquier, and S. Plaza, in [Tj, modified it to a 
third order approximation function Mf(x ) = Nf(x) — f (Nf(x)) / f (x) that is free of sec¬ 
ond derivatives and shows a remarkable robustness when compared to other methods. 
On the other hand, it is known that the classical Newton’s Approximation Function 
Nf, when considered as a discrete dynamical system, shows chaotic behaviour, at least 
from two possible acceptions of the concept of chaos: first in the sense of T-Y. Li and 
J.A. Yorke [5], that is, existence of periodic points of any prime period, and second 
in the sense of R. Bowen [3], that is, a strictly positive topological entropy, as shown 
by M. Hurley and C. Martin [3], see also D.G. Saari and J.B. Urenko [7] for similar 
investigations. 

In p], the chaotic behaviour of Mj was numerically pointed out for polynomials 
of order less or equal than 3 by using a bifurcation diagram similar to that of the 
logistic map. They left open the question of performing a qualitative analysis on 
the discrete dynamical system associated to Mf in order to mathematically prove its 
chaotic behaviour. The main result of this article is Theorem 13.31 that shows that for 
functions / of Newton type, see Section [3] for definition, with at least four roots and 
infinite limits of opposite signs at ±oo, Mf has periodic points of any prime period and 
the set of points a at which the approximation sequence (M"(a )) ne n does not converge 
is uncountable. For example, when considering polynomials, this result applies to all 
odd degree polynomials with a certain distribution of real roots. In view of Lemma [2131 
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Theorem 13.31 might be extended to other classes of Newton’s functions / having at least 
four roots and for which Mf has at least two bands that cover the whole real line R, 
that is, there are two pairs of disjoint open intervals, formed by consecutive critical 
points of /, that are mapped by Mf onto the whole R, see [0] for precise terminology. 

In addition, in Theorem 14.11 we observe that in the Scaling Theorem from [T], which 
says that the dynamics of Mf is stable under affine conjugacy, as well as in its damped 
version as in |2], analyticity of / can be replaced by its differentiability. The Scaling 
Theorem is essential for the analysis performed in [1| because it reduces the study of 
the dynamics of Mf for a general class of functions / to the study of the dynamics 
of Mf for a considerably smaller class of simpler functions. For example, in order to 
understand the dynamics of Mf for quadratic polynomials /, it suffices to study only 
the dynamics of Mf for the quadratic polynomials x 2 , x 2 + 1 , x 2 — 1 , while for cubic 
polynomials /, it suffices to study only the dynamics of Mf for the cubic polynomials 
x 3 , x 3 + 1, x 3 — 1, x 3 + yx + 1, with 7 e R. 

We might also add results on the lower estimation of the Bowen’s topological entropy 
for Mf but, once Theorem 13.31 is obtained, these follow in an almost identical fashion 
as in [4j, when considered the two kind of bands induced by the critical points of / 
for Mf. Also, we observe that, for the class of functions that make the assumptions 
of Theorem 13.31 the damping with a parameter A, cf. [2], does not change either the 
chaotic behaviour expressed as the existence of periodic points of any prime period or 
the uncountability of the set of all real numbers a for which the iterative sequences 
(M\, f( a ))n=i diverges. 


2. Preliminaries 

In this section we collect results related to periodic points of continuous real func¬ 
tions. We start by proving a sequence of lemmas that are essentially contained in [4J 
and more or less implicit/explicit in the works of N.A. Sharkovsky [ 8 ] and T.-Y. Li and 
J.A. Yorke pj. 

The first lemma is a well known fixed point result and a direct consequence of the 
Intermediate Value Property for continuous functions. 

Lemma 2.1. If I and J are compact intervals and f: I J is a continuous function 
with /(/) D I, then f has a fixed point. 

The next lemma is important for understanding the dynamics of continuous real 
functions, e.g. see Lemma 0 in [5]. We provide a proof for consistency. 

Lemma 2.2. Let J and K be nonempty compact intervals and let f: J —» R be a 
continuous function such that f(J) 5 K. Then, there exists a nonempty compact 
interval L C J such that f(L) = K. 

Proof. Let K — [a, b\. If a = b then we apply the Intermediate Value Theorem and 
get c G J such that /(c) = a and let L = [c, c], so let us assume that a < b. Since 
the set {x E J \ f(x) = a} is compact and nonempty, there is a greatest element c 
in this set. If f(x) = b for some x > c with x C J, then x > c and letting d be the 
least of them, by the Intermediate Value Theorem f([c,d\) = K and we let L = [c,d\. 
Otherwise, f(x) = b for some x < c and let c! be the largest of them. Then let d' be 
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the smallest of the set of all x > d with f(x) = a so that [c',d'\ is an interval in J. 
Then f([cf, d']) = K and we let L — [c', d']. □ 

The main technical fact we use is a refinement of Lemma 2.2 in [4j, implicit in the 
proof of Sharkovsky’s Theorem [8], see also [5|. Recall that, a point a G M is called a 
periodic point for a function g: M —> M if there exists n G N such that g n (a) = a, and 
the least n G N with this property is called its prime period. 

Lemma 2.3. Let g: M. —> R be a function and let Ii,l 2 , ■ ■ ■ ,Ik be compact, disjoint 
and nondegenerate intervals, with k > 2, such that, for all m G {1,2 g is 
continuous on I m and 

k 

g(Im) If [_J Ij- 

3 =1 

Then: 

(a) For each n G N, g has at least k(k — l) n_1 periodic points of prime period n, in 
particular, g has periodic points of any prime period. 

(b) The set of all real numbers a for which the orbit (g n (a)) nG fq makes a sequence 
that does not converge is uncountable. 

Proof, (a) For any n G N, take any sequence of length n with j l G {1, 2,..., k) 

and let j n+ 1 = j\. Considering the sequence of intervals by assumption, 

k n 

5(4) 2 U JJ-2U4 2 4- 

j=1 i=l 

By Lemma 12.21 there exists a compact interval A n C I n such that g{A n ) = Ij 2 . If 
n = 1, we observe that, since g(A ]l ) = I ]2 = I n , by Lemma [2.II it follows that g has a 
fixed point in the compact interval A n . 

If 77, > 2, then 

k Tl 

9 2 (A.) = 9(4) 2 U Jj 2 U4 2 4. 

7=1 *=1 

and by Lemma 12.21 applied to (? 2 we obtain a compact interval A,- 2 C A n C I n such 
that g(Aj 2 ) = Ij 3 . Proceeding in a similar fashion, after n consecutive applications of 
Lemma [2.21 we obtain compact intervals 

Aj cl Aj . c ••• c Aj, cl I j 

such that 

(2.1) g\A,f) = 4 +i , i = l,...,n, 

in particular, g n (Aj n ) = / >n+1 = I n . Then, by Lemma 12 . II there is a fixed point a G A, n 
for g n , hence a is a periodic point for g of period n. 

We observe now that, for n — 1 there are exactly k different choices for j\ — 1,..., k 
and, since the intervals Ji,..., Jk are mutually disjoint, the k fixed points obtained, as 
explained before, are all different, hence g has at least k fixed points. 

For n > 2, if j\ ^ ji for all i — 2,..., n, then the periodic point a for g obtained before 
from the sequence (ji)” =1 has prime period n. Indeed, for all 2 < i < n, a G Aj. C A Jz _ 1 
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and, by (12.11) it follows that g l ~ l {a) G Ij i hence, since I n fl Ij i = 0 it follows that a 
cannot have any period less than n. 

We consider now two different sequences (j,;)” =1 , {k)2=i as before, hence there exists 
r G {1, 2,..., n) such that j r j- l r . If the hxed point a for the two sequences of intervals 
(Iji)i =l an d (4)?=i, obtained as before, is the same, then (g l (a))™ =0 , the orbit of a under 
g, is the same for both sequences. But then, g r ~ 1 (a) G Ij r D Ii r — 0 hence we have 
a contradiction. Therefore, for the k(k — l )" -1 different sequences (h)” =1 of length n, 
formed with elements from the set { 1 ,..., k} and subject to the condition j\ ^ ji for 
all i = 2,... ,7i, we have k{k — l) n different hxed points of prime period n. 

(b) For any infinite sequence (ji)°Z 1 with elements from {1, 2,..., k}, which is not 
eventually constant, by the same construction as above we get a sequence of nonempty 
compact intervals (db,'. ( ) 4 i subject to the properties 

(2.2) • • • A Ji+1 C A jt C • • • C A ja C A h C I n , 
and 

(2.3) g ^ Aj .) = 4 +1 , iGl 
By the Finite Intersection Property we have 

OO 

A =f] A i< t 6 

2—1 

hence, any point a G A has the property that its orbit (p*(a ))“ 1 makes a sequence that 
does not converge. The sequence (^(a))^ does not converge since, for all i G N we 
have g l (a) G I M+l , the sequence with elements from {1, 2,..., k} is not eventually 

constant, and the compact intervals ii,..., /*, are mutually disjoint. 

Let us consider two different sequences and (/j)“ 1 , formed with elements from 

the set { 1 ,..., k} and not eventually constant. As before, to the sequence (jj )“ 1 we 
associate the sequence of nonempty compact intervals (A,-.)°f x subject to the properties 
(12.21) and (12.31) . and let a G f)4i 4/V Similarly, there is a sequence (Bi i ) c *L 1 of nonempty 
compact intervals subject to the properties 

■■■B li+1 CB h C...CB l2 CB h CI jl , 

and 

g\B h ) = i kw ■€*, 

and let b G H4i B^ ^ 0. We claim that a ^ b. Indeed, there exists r G N such that 
j r 7 ^ l r , hence g r (a ) G Ij r D Jp = 0, a contradiction. 

In conclusion, there are as many real numbers a for which the sequence (g l (a))jZi does 
not converge at least as many as sequences (ji)jZi formed with elements from the set 
{ 1 ,..., k} and that are not eventually constant, and the latter set is uncountable. □ 

3. The Dynamics of M f 

Following [0], the Newton Class is defined as the collection of all real functions / 
subject to the following conditions: 

(nfl) / is of class C 2 (M). 

(nf2) If f(x ) = 0 then f'(x ) ^ 0. 

(nf3) If f(x) = 0 then f"{x) ^ Q. 
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Remarks 3.1. Let / be a Newton map. 

(a) Clearly, both / and its derivative f do not have multiple roots. 

(b) The roots of /, and similarly the roots of /', do not have finite accumulation 
points. Indeed, if (x n ) ne ^ is a sequence of distinct roots of / that accumulates to some 
io 6 I then, by the Interlacing Property of the roots of / and its derivative /', it 
follows that there exists a sequence (c n ) ne N of distinct roots of f converging to xo- 
Since both / and f are continuous, it follows that xq is a root for both / and /', 
contradiction with the statement at item (a). A similar argument shows that the roots 
of f do not have finite accumulation points. 


For a differentiable function / on an open set PCI, the Classical Newton’s Ap¬ 
proximation Function Nf is the function defined for all x G D such that f\x ) 7 ^ 0 
by 


(3.1) 


N f (x) 


_ /M 

/'W 


Following p], the Modified Newton’s Approximation Function Mf is the function, de¬ 
fined for all x G D such that f'(x ) 7 ^ 0 and Nf(x) G D , 


(3.2) 


M f {x) 


f(x) 

f(x) 


f(x 


f 0 ) 1 
F(*) J 


f'(x) 


or, in terms of Nf, 


(3.3) 


M,(x) = N,(x) 


WM) 

f(x) 


The following lemma provides some information on the behavior of Nf in the neigh¬ 
bourhood of the critical points of a Newton map /, see Remark 1.3 in [4]. 


Lemma 3.2. Let f be a Newton map and let C\ and c 2 be two consecutive roots of f 
such that in the interval ( 01 , 02 ) there is a unique root of f. Then 


lim Nf(x ) = — lim Nf(x) = ± 00 , 

X— >ci+ X—>C2 — 


Proof. Indeed, near c\+ and c 2 —, / has opposite signs since it has a unique root inside 
the interval (ci, c 2 ), while f does not change its sign and goes to zero when approaching 
both ciT and c 2 —. Since the term x is majorised by the term near ci+ and c 2 —, 
we have the result. □ 


Here is the main result that shows that the Modified Newton’s Approximation Func¬ 
tion Mf provides chaotic behaviour for a large class of Newton’s functions /. 

Theorem 3.3. Let the function f: M — )• M have the following properties: 

(i) f is a Newton’s function. 

(ii) lim f(x) = — lim f(x) = ± 00 . 

ir—>-+00 x — y —00 

(iii) / has at least four real roots. 

Then: 

(a) Mf has periodic points of any prime period. 
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(b) The set of all real numbers a for which the sequence (Mf (a)) nG N does not converge 
is uncountable. 

Proof. By Rollc’s Theorem, between any two consecutive roots of / there exists a 
root of its derivative f hence, by remarks 13.11 and property (iii), we can choose four 
consecutive roots rq < r 2 < r 3 < r 4 and three or four roots C\ < C 2 < c ' 2 < C 3 of f in 
the intervals (ri,r i+ i) for i = 1,2,3, respectively, with a unique root of / inside each 
of the intervals (ci, C 2 ), {c 2 , c 3 ) and with no other root of /' there. With this choice, f 
does not change its sign on the intervals (ci, C 2 ) and (c 2 , c 3 ). By using Lemma [3721 and 
property (ii), we have 

(3.4) 

In the following we show that 


lim f(Nr(x )) = — lim f(Nr(x )) = ± 00 . 

rr—>-ci+ x^C2 — 


(3.5) 


lim Mf(x ) = — lim Mf(x ) = ± 00 . 

X— »C1+ X—>C2 — 


Indeed, when x —* c 3 +, the first equality will follow if we show that 

lim = -_m zmm = ±COi 

z->ci+ f'[x) 

and, taking into account that lim f{x) = f(cf) G R, we observe that the latter will 

X—>C\ + 

follow if we prove that 

-f(NAx )) 

lim , - 777—\ - = ±CX) ’ 

Hii+ J (X) 

which actually follows from (j3.4[i and the fact that f does not change its sign in the in¬ 
terval (ci, C 2 ). Hence, the first equality in (13.51) is proven. The fact that lim Mf(x ) = 

X—>C 2 — 

±00 is proven similarly, the only thing that remains to be shown is that the two limits 
in (13.51) have different signs, which actually is a consequence of (j3.4[) and the fact that 
f has constant sign on ( 01 , 02 ). Therefore, (13.51) is proven. 

Similarly, we have that 

lim MAx ) = — lim Mt{x ) = ± 00 . 

x—tc' 2 + x ^c 3 

Therefore, there exists e sufficiently small such that, letting I\ = [c\ + e, c 2 — e] and 
I 2 = [c' 2 + e, c 3 — e] we observe that I\ D I 2 — 0 and that Mf(Ij) D [ci, c 3 ], for j = 1 , 2 , 
hence Mf{If) and contain I\ U 12 - Finally, Lemma [2.31 is now applicable with 

k — 2 , which finishes the proof. □ 

Remark 3.4. In order to reduce the chaotic behaviour and improve numerical param¬ 
eters of approximation for lower order polynomials, in [2] and [ 6 ] damped Newton’s 
methods have been considered. More precisely, letting A be the damping parameter, 
one defines N X j and M\j as follows: 

f(x) 


(3.6) 
and 

(3.7) 


N\, f (x) = x - A 


= N X j(x) - A 


f'{x) ’ 

f(N XJ (x)) 

f'(x) 
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It is easy to observe, by inspection, that Lemma 13.21 and Theorem 13.31 remain true 
if Nxj and replace Nf and, respectively, My, for arbitrary damping parameter 

A > 0, hence the chaotic behaviour characterised by existence of periodic points of 
any prime period, as well as the uncountability of the set of points of divergence of 
iteration of My, remain unaltered by damping, for the class of functions considered in 
Theorem 13.31 


4. The Scaling Theorem 


In this section we observe that the Scaling Theorem, cf. Theorem 1 in |lj, which says 
that My is stable under affine conjugation, remains true if we replace the analyticity 
condition on the function / is with its differentiability. 

Theorem 4.1 (The Scaling Theorem). Let f be a differentiable function on M. Let 
T(x) — ax + b be an affine map with a f 0. Then, for all x G M with f'fx ) f 0, we 
have 

(T o M foT o T” 1 )(x) = M f (x). 

We first prove a lemma saying that the transformation Nf satisfies a similar property 
of stability under linear conjugation. 


Lemma 4.2. Under the assumptions as in Theorem f.L for all x el with f'(x) 0, 

we have 

(T o N foT o T- 1 ) (x) = N f {x). 

Proof. Clearly, for all y e M we have (/ o T)'{y) = af'{ay + b ) hence, if f'[ay + b) 0 
we have 

/ {ay + b ) 


(4.1) 

and then 


N f°T{y) = y 


af(ay + b )’ 


— N f (ay + b) = (N f o T){y). 


Since T x (t) = —-, letting x = ay + b we get 


(T o Nj oT o T~ 1 )(x) — JVf(x). 


□ 


Proof of Theorem f.l. By fj3.2jl . for arbitrary j/gR with f'{ay + b) f 0, from (14.11) we 
obtain 

(T O MfoT) (y) = aNf oT (y ) + b - ^ aN f oT ^ + ^ 


and then, letting x = ay + b, 

(:T O M foT O T- 1 ) (x) = (T O N foT O T" 1 ) (x) - 


f'{ay + b ) 

/((To N foT o T -1 ) (x)) 
f'(a(^) + b) 
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whence, applying Lemma 14.21 

= N,(x) - hMld = Mf (x). □ 

Remark 4.3. The Scaling Theorem remains true, with almost exactly the same proof, 
for the damped Newton’s function My/, see (13.71) . for arbitrary damping parameter 
A 7 ^ 0, more precisely, we have 

(ToM XJoT oT- 1 )(x) = M xj (x), 

for all i 6 R such that f'{x) ^ 0. For the case of the damped Newton’s function N\j, 
see (13.61) . the corresponding generalisation of Lemma [4.21 

(ToN XJoT oT- 1 )(x) = N xj (x), 

follows from the proof of Theorem 2.1 in pa]; although it was stated for analytic func¬ 
tions /, the assumption was not used in that proof. 
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